Superfluid phases of the three-species fermion gas 
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We discuss the zero temperature phase diagram of a dilute gas with three fermionic species. We 
make use of solvable limits to conjecture the behavior of the system in the "unitary" regions. The 
physics of the Thomas-Efimov effect plays a role in these considerations. We find a rich phase 
diagram with superfluid, gapless superfluid and inhomogeneous phases with different symmetry 
breaking patterns. We then discuss one particular possible experimental implementation in a system 
of ^Li atoms and the possible phases arising in this system as an external magnetic field is varied 
across three overlaping Feshbach resonances. We also suggest how to experimentally distinguish the 
different phases. 



Cold atomic gases provides are great laboratory for 
many-body physics. This is due to the variety of atomic 
isotopes that can be trapped and the tunability of the 
atomic interactions. Another point of interest is the "uni- 
versality" of some of their properties: since the detailed 
shape of the interatomic potential is not probed in dilute 
systems, many different kinds of interactions can be ap- 
proximated by a contact potential and different systems 
present analogous characteristics. 

Our goal in this letter is to explore the phase dia- 
gram of a dilute ultracold gas composed of three fermionic 
species [1]. It can be considered a generalization of the 
well studied two fermionic species gas and brings about 
some new features like competition among interactions 
and the consequences of the Efimov efect [2] Let us first 
summarize the main properties of the two-fermion case 
as they will be useful in the more involved three-fermion 
gas. There are three relevant distance scales in this sys- 
tem: ro measures the range of the interatomic forces, 
a is the s-wave scattering length and the typical 

interatomic distance (n is the particle density). We will 
consider only the cases where rg <C a, 1/ n^/^ since only in 
this case the results have some degree of universality. The 
cases where n\af <gC 1, where an expansion in powers of 
n\a\^ is possible, are more easily analyzed. In the regime 
na^ <^ 1, a > the atoms bind into diatomic molecules 
much smaller than the separation between them while in- 
teracting weakly with each other. These molecules then 
condense to form a weakly coupled superfluid. If a < 0, 
though, there are no molecules. Still, the atoms form 
Cooper pairs and the ground state of the system is a BCS 
superfluid. In the intermediate regime — 1 ^ 1/na? ^ 1, 
the so-called "unitary" region, the absence of a small pa- 
rameter complicates the theoretical analysis but one can 
still infer some of its qualitative properties. Notice that 
this region is scaled down (in energy) version of the di- 
lute neutron gas, exemplifying the universality pointed 
out above. In both the BEC limit 1/na? ^ 1 and the 
BCS limit \/na? ^ — 1 the condensation of molecules or 



Cooper pairs breaks the wave function phase symmetry. 
That is, if we denote by "01 ("02), the annihilation op- 
erator for the atom of type 1 (2), the symmetry of the 
hamiltionian ^l)i^2 e*"V'i,2 is not a symmetry of the 
ground state due to the non- vanishing value of (V'iV'2)- 
Since the symmetry breaking pattern is the same on both 
sides of the "unitary" region there is not reason to expect 
a qualitative change (a phase transition) in between the 
BEC and the BCS limits. This picture is vindicated both 
by experimental [3] and numerical [4] results. 

Similar arguments can be made in the three-fermion 
case of interest here. There are two main new ingre- 
dients absent in the two-fermion case. The first is the 
competition between the pairing/binding among different 
atom pairs. The second is related to three-body corre- 
lations. Ordinarily, two-body correlations dominate the 
physics of dilute systems but that may not be true in 
the regime rg/a <C 1. Indeed the Thomas effect [5] (the 
collapse of the three-body system as ro ^ at fixed 
a) and the Efimov effect [2] (the accumulation of three- 
body bound states at threshold as a — > oo at fixed rg) 
are consequences of the strong attraction among three 
resonating particles. These two closely related phenom- 
ena [6] can occur only if the three particles in the system 
can overlap in space, which is forbidden in the case of 
only two fermion species. The sensitivity to the physics 
of three overlapping particles brings out a dependence of 
the low energy processes in a new parameter [2, 7] un- 
related to two-body scattering which, in the language of 
Effective Field Theories, is the value of a contact three- 
body force [8]. Currently, the experimental observation 
of the Efimov effect is a subject of great interest due to 
the fundamental aspects it represents for the three-body 
physics in ultracold quantum gases. Recently an analy- 
sis of the three-body recombination [9] in Bose-Einstein 
condensates have shown evidence of this intriguing effect 
[10]. 

At zero temperature, with equal densities of the three 
species, the phase diagram of the system is a function 
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FIG. 1: Conjectured phase diagram. The horizontal bands 
(blue online) show the narrow areas of atom-molecule sepa- 
ration. The solid (red online) represent the phase boundary 
between the Z2 and SU{2) x f/(l)Z2 phases. The dashed line 
is where a Fermi surface appears (Lifshitz transition). 



of four dimensionless variables. The dependence on the 
three-body parameter is hmited so we keep it fixed. The 
remaining three parameters can be taken to be na'fj , {aij 
is the scattering length between particles i and j and 
n the density of each atomic species). We will initially 
restrict our discussion to the slightly idealized situation 
where two of the scattering lengths are equal, thus as- 
suming ai2 = = a, 031 = a. Our strategy to map the 
zero temperature phase diagram as a function of a and d 
will be to understand the system in the soluble limits of 
the phase diagram where n|ap,n|ap <g; 1 and make ed- 
ucated guesses about how to interpolate between these 
limits. This is the same successful strategy discussed 
above for the two-fermion case. 

The system can be described by the hamiltonian 
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+ ff'(IV'i|>2|' + |V'2n^3|') + GIV^iH^anV-al'. 

The couplings g, g and G depend on a cutoff and are 
determined respectively by a, a and, in the case of G, 
on a three-body observable like one of the atom-dimer 
scattering lengths for some value of a, d. The hamiltonian 
in eq. (1) has a SUi3{2) x J72(l) x t/i23(l) symmetry, 
where the first factor acts only on atoms of type "1" and 
"3" , the second on atoms of type "2" only and the last 
is a common phase change in all atom fields. 

Let us first consider the case where both scattering 
lengths are positive and na'^^nd^ ^ 1, located in the 
third quadrant in Fig. 1. There all three pairs can form 
bound states of energy B = 1/Ma? or B = 1/McP and 
whose spatial extent is much smaller than the interpar- 
ticle distance. If we denote the number of bound states 
between particles i and j by riij and the number of un- 
bound atoms of type i by rii we have that, to leading 
order in na^ and no'^, the energy of the system is given 
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which is just the sum of the molecule binding energies and 
the energy of the Fermi gas of unbound atoms. We dis- 
regarded the interaction between unbound atom-dimer, 
unbound atom-unbound atom and dimer-dimer as these 
contribution are, as long as their respective scattering 
lengths are of the order of ~ a, a, suppressed by powers of 
na^,nd'^ <C 1 . As we will see later the atom-dimer scat- 
tering length may diverge for some very particular values 
of the parameters as a consequence of the three-body dy- 
namics and will be necessary to amend this analysis. We 
find the ground state by minimizing eq. (2) with the re- 
striction that the total number of atoms of each species is 
n (ni -|-ni2 +n3i = n and similarly for the other species). 
We find three different cases: 
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This result is physically clear. In the first case, it is favor- 
able for all atoms to form bound states, and n/2 dimers 
are formed between each species pair. All dimers species 
condense leading to a non-zero value of the ground state 
expectation value (ipiijjj) for all This condensate 

breaks the symmetry from SUi3{2) x ?72(1) x ?7i23(l) 
down to I12, corresponding to a common sign flip on 
all atom fields. If B increases (i.e., d decreases) be- 
yond B > 2B it is favorable to break some of the "12" 
and "23" dimers in order to increase the number of "31" 
dimers, leaving some "2" atoms unbound. The conden- 
sates (■01 "02) and (0203) will decrease in value but the 
symmetry breaking pattern will remain the same. The 
formation of a Fermi surface for the "2" atoms indi- 
cates a topological (Lifshitz) phase transition (denoted 
by the dashed line in fig. (1)). As we keep increas- 
ing B (decreasing d) there comes a point where all "2" 
atoms are unbound and only "13" dimers can be found. 
At this point (V'i'02) = (0203) = and the symmetry 
breaking pattern becomes SUi3{2) x U2{1) x ?7i23(l) 
SUi3{2) X U2{1) X Z2. The change in symmetry of the 
ground state indicates the presence of a phase transition, 
marked by the full line in the third quadrant of fig. (1). 

We now consider the opposite side of the phase dia- 
gram, i.e., the region where a, a < 0, and n|ap, n\d\^ ^ 1 
(first quadrant in fig. 1). In this region there are no two- 
body bound states, however, symmetries can be broken 
by the condensation of Cooper pairs. The most general 
s-wave condensate can be written as 
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where A = (A23, A31, A12) is a complex vector describ- 
ing the fermionic gaps. In the regime considered the BCS 
theory (mean field) can be used. Using the methods of 
reference [12] wc find the thermodynamical potential in 
this approximation to be 

1 r rl^k , 



IA12P + IA23P 
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where €{k) = k^/2M — ji and /it is the chemical poten- 
tial. In cq. (4) wc traded the couplings constants by 
the scattering lengths. In general, in order to maintain 
equal densities for the three species, different values for 
the chemical potentials for the three species arc required. 
At low density though, those differences (of the order of 
A^/kp, where A is the gap and kp the Fermi momen- 
tum) are suppressed by the exponentially small size of 
the gap and will be neglected here. The ground state is 
found by minimizing eq. (4). The first term in eq. (4) 
is a function only of the combination |Ap which is a 
'S'C^i23(3) invariant. The remaining terms don't have this 
symmetry and, consequently, tip the balance in favor of 
one direction in A'^ space. The location of the minimum 
then depends on the competition between the strength 
of the interaction between the atoms, that is, on relative 
sizes of a and a. If \a\ > \a\, the energy is minimized by 
having A13 7^ and A12 = A23 = 0. The atoms of type 
"1" and "3" pair up with each other while the atoms of 
type "2" remain impaired. The symmetry breaking pat- 
tern is 5t/i3(2)xC/2(l)xi7i23(l) ^ S'C/i3(2)xC/2(l) xZ2. 

In the opposite case, \a\ < \a\, the energetically fa- 
vored state has A13 = and |Ai2p + IA23P 7^ 0. It 
is equally favorable for A12 or A23, or a combination 
of them to be non-zero, due to the SUi3{2) symmetry 
of the hamiltonian. The symmetry breaking pattern is 
SUi3{2) X C/2(l) X [/i23(l) — > Z2. In this phase atoms 
of type "1" and "3" do not pair up with each other but 
pair up with atoms of type "2" . As a consequence of the 
different pairing patterns we conclude that there must 
be a phase transition along the line a = a in the first 
quadrant. Both phases are superfluid since both contain 
Cooper pairs but the 5J7i3(2) x J72(l) x Z2 phase also has 
metal-like properties on account of the unpaired atoms 
of type "2". 

Let us consider now the second quadrant in (Fig. 1), 
approaching it from the third quadrant by moving along 
the vertical direction. As we saw, in the third quadrant 
just below the horizontal axis all atoms are bound. As we 
cross the horizontal axis, atoms of type "1" and "3" are 
no longer bound but can still form Cooper pairs. In fact, 
the transition between bound atoms and Cooper pairs 
does not change the symmetries of the ground state (that 
remains Z2) and no phase transition is expected. This is 
perfectly analogous to the BEC/BCS transition discussed 



of the two-fcrmion system. Finally, let us consider the 
fourth quadrant. There is also a smooth connection be- 
tween third and fourth quadrants as a flips sign. In the 
third quadrant, just to the left of the vertical axis, atoms 
"2" are attracted to the other ones, but all atoms "1" and 
"3" are already bound and not available for binding with 
"2" and the attraction to atoms "2" is not strong enough 
to force a rearrangement of the binding. As we move 
horizontally towards the right and cross the vertical axis, 
the interaction of atoms of type "2" with atoms "1" and 
"3" becomes weaker and does not even support a bound 
state. Consequently, no qualitative change occurs and 
no phase transition is expected. The symmetry break- 
ing pattern remains the same and no phase transition is 
expected. 

Up to now we have taken into account only two-body 
correlations - the existence of two-body bound states 
and/or Cooper pairs - to discuss the ground state of 
the system. This is reasonable in dilute systems with 
short range interactions as the probability of having 
more particles at the same point in space is small. Our 
previous discussion receives corrections proportional to 
n|a|'^,n|ap <C 1 and is, consequently, robust. This ar- 
gument however, implicitly assumes that parameters de- 
scribing the three-body physics, for instance, the scat- 
tering length between an atom and a molecule, are of 
the same order as a or a. The presence of a three-body 
scale , however, may invalidate the argument. This is in- 
deed what happens for some specific values of the pa- 
rameters. It is known that the atom-dimer scattering 
length becomes arbitrarily large for specific values of the 
parameters. When the atom-dimer interaction is strong 
and repulsive the possibility of spatial separation between 
atoms and molecules arises. Free atoms and molecules 
appear simultaneously only in the third quadrant of the 
phase diagram and only the atom of type "2" is found 
unbound among molecules of typ(^ "13" . The criterion for 
the (linear) stability against atom/molecule separation in 
a boson-fermion (atom-molecule) mixture is [13] 
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where ai3_i3 Ri 0.6a [14] is the scattering length between 
two "13" dimers and, for a < a we have [2] 

ai3,2 = a{A + Btan{solog-^ + ^)), (6) 

where A, B and so » 1.006 are numerical constants. Un- 
fortunately, $ is a function of the three-body force G 
and thus can only be determined with the knowledge of 
a three-body observable. For a < a, there is no enhance- 
ment and ai3.2 ^ a. Notice that a such divergence of 
the atom-dimer scattering length (the poles in Eq. (6)) 
only occur for fermionic gases if there is, at least, three 
different species. In other cases the Pauli principle for- 
bids the atoms to overlap and limits the strength of the 
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interaction among the three particles. For a near one 
of the poles of Eq. (6) the stability condition given by 
Eq. (5) is violated and atoms of type "2" separate from 
the "13" molecules. The experimental observation of this 
phenomenon, as well as its repetition for different poles in 
^13, 2, would provide a striking confirmation of the Efimov 
effect. The bands where this instability occurs are re- 
peated every time a is increased by a factor e^/^° 22.7. 
We have indicated two of these bands in our phase dia- 
gram in Fig. 1 . Their fate as we cross the vertical axis is 
unclear, i.e., they may or may not extend into the fourth 
quadrant. The precise position where these bands occur 
when we approach the unitary regime {—l/n\a\^ <C 1) 
varies for each systems as it depends on three-body in- 
formation encapsulated in G. It is hard to determine 
whether this repetition continues to very small values of 
— l/n|a|^ and we then enter the unitary regime where our 
methods fail. 

Our considerations can be easily extended for the case 
where all three atom-atom scattering lengths are differ- 
ent. Here we will consider the phases arising in a con- 
crete possible realization of the three- fermion gas. It 
was predicted in Ref. [11] the existence of three Feshbach 
resonances between atoms in the lowest three hyperfine 
states of ^Li for values of the magnetic field B in the 
range of 60 mT < B < 100 mT. In this case, all scat- 
tering lengths are positive for small fields and negative 
for large fields. The resonances occur around B « 69 
mT ("13"), B ?» Sl.lmT and ("12") and B « 83.4 mT 
("23"), where the respective scattering length diverges 
and flip sign. For the largest B in this range all 
are negative, the ground state of the system is an ana- 
logue of the phase in the first quadrant in fig. (1) and the 
condensate is driven by the largest scattering length in 
absolute value. A (V'iV'2) condensate forms which breaks 
the Ui{l) X 1/2(1) X ^73(1) symmetry of the hamiltonian 
down to t/i_2(l) X Us{l) X Z2 (?7i_2(l) corresponds to an 
equal but opposite phase rotation of the atoms "1" and 
"2"). As the field is lowered below B = 81.1 mT, 012 
and a23 flip sign but that initially does not change the 
phase of the system. Atoms "1" and "2" that formed a 
Cooper pair at large B now form a bound state. Lower- 
ing B further still, at some point the binding energies for 
the pairs "12" and "23" are close to each other and large 
enough to dominate the energetics of the system. At this 
point, type "12" and "23" molecules form and condense. 
The creation of condensates (V'iV'2) and (^2^3) break the 
symmetry down to Ui{l) x 1/2(1) x 1/3(1) — »■ Z2, indicat- 
ing the existence of a different phase from the one at high 
values of B. This is the analogue of the phase transition 
in the third quadrant of fig. (1). As we mentioned above, 
this picture may be modified by the atom-molecule sep- 
aration for certain values of B below the "12" resonance, 
but we cannot predict for which value of B this will oc- 
cur due to our ignorance of the parameter $ (or G) . The 
phase separation would be observed in practice if one 



would be able to scan ai2 through a factor of e'^/®" « 22.7 
in order to be sure to cross one of instability bands. A 
measurement of three-body observables involving these 
three states of ^Li would determine G (and $) and would 
allow us to make a more quantitative prediction. 

We close by pointing out that the two different su- 
perfiuid phases can be distinguished by the properties of 
their vortices. In the three ®Li system, for instance, and 
at high fields, atoms of type "3" are unpaired and should 
not participate in in the formation of superfiuid vortices. 
For lower fields, below the point where the phase transi- 
tion occur all atom types can only rotate by the formation 
of vortices. The creation of these vortices by means of 
probes coupling to only one atom type is a direct way of 
exploring the phase diagram we have discussed. 
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